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Abstract. We examine here the classical dynamics of cold atoms in square optical lattices, 
i.e. lattices obtained with two orthogonal stationary plane waves. Contrary to much of the past 
studies in this domain, the potential is here time independent and non dissipative. We show that, 
as a function of the experimental parameters, very different behaviors are obtained, both for the 
dynamics of atoms trapped inside individual sites, and for atoms travelling between sites: inside 
the sites, chaos may be a main regime or, on the contrary, may be negligible; outside the sites, 
chaos sometimes coexists with other regimes. We discuss what are the consequences of these 
differences on the macroscopic behavior of the atoms in the lattice, and we propose experimental 
measurements able to characterize these dynamics and to distinguish between the different cases. 

1. Introduction 

The cooling of atoms to extremely low temperatures, through the use of magneto- 
optical traps (MOT) , opened since the mid eighties fantastic possibilities to increase 
our experimental knowledge of the quantum world. The most spectacular realization 
was the achievement of the Bose-Einstein condensation, and thus of macroscopic 
quantum objects. However, even in the classical world, the possibility to study the 
dynamics of atoms not "blurred" by the Doppler effect is very exciting. This requires 
to develop tools to manipulate the atoms, for e.g. guiding them or "designing" their 
phase space. 

Optical lattices provide such tools: their versatility allows to manipulate atoms 
with an extreme precision and a relative ease (Guidoni & Verkerk 1999). Because 
of these qualities, they represent an outstanding toy model, and have recently 
attracted increasing interest in various domains. Condensed matter systems and 
strongly correlated cold atoms in optical lattices offer deep similarities, as in the 
superfluid-Mott insulator quantum phase transition (Greiner et al. 2002), or in the 
Tonks-Girardeau regime (Paredes et al. 2004). Here, the interactions between atoms 
play a crucial role, and require the use of a Bose-Einstein condensate. In particular, 
instabilities are expected in the Gross-Pitaevskii equation, because of the non-linear 
term (Thommcn et al. 2003, Fang & Hai 2005, Kuan et al. 2007). Quantum com- 
puting also requires a coupling between atoms: the optical lattices appear to be 
an efficient implementation of a Feynman's universal quantum simulator (Jaksch & 
Zoller 2005), and are among the most promising candidates for the realization of a 
quantum computer (Mandel et al. 2003, Vollbrecht et al. 2004). On the other hand, 
noninteracting atoms also exhibit interesting behaviors. In this case, the physics is 
essentially that of a single atom. A higher number of atoms simply increases the 
observable signal. That is the case in statistical physics, where cold atoms in optical 
lattices, through their tunability, made possible the observation of the transition 
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between Gaussian and power-law tail distributions, in particular the Tsallis distri- 
butions (Douglas et al. 2006, Jersblad et al. 2004). or that of Anderson localization 
(Billy et al. 2008, Roati et al. 2008, Chabe et al. 2009). 

Non interacting cold atoms appear also to be an ideal model system to study the 
dynamics of a system in its classical and quantum limits. Both are closely related, 
as the latter is only defined as a function of the former. In particular, quantum 
chaos is defined as the quantum regime of a system whose classical dynamics is 
chaotic. A good understanding of the classical dynamics is therefore an essential 
prerequisite to the study of quantum dynamics. In non dissipative optical lattices, 
both the classical and the quantum limits are experimentally accessible, and it is 
even possible to change quasi continuously from a regime to the other (Steck et al. 
2000). Moreover, the extreme flexibility of the optical lattices makes it possible to 
imagine a practically infinite number of configurations by varying the complexity 
of the lattice and the degree of coupling between the atoms and the lattice. 

Many results have been obtained during these last years in the field of quan- 
tum chaos (Steck et al. 2000, Lignicr et al. 2005). However, all these works have 
used very simple potentials, mainly ID. Chaos is obtained only with a periodic 
(or quasi-periodic) temporal forcing of the depth of the lattice (Steck et al. 2000, 
Lignier et al. 2005), and only the temporal dynamics of the individual atoms is 
studied. The introduction of this external clock and the restriction to ID potentials 
reduce considerably the generality of these results and the type of possible dynam- 
ics. In particular, the behaviors related to the appearance of new frequencies or 
to a frequency shift (quasi-periodic and homoclinic bifurcations, for example) are 
impossible. 

If we want to break these limitations, several problems have to be examined: 
what type of time-independent lattice will lead to a reasonably complex dynamics? 
What are the relevant quantities to characterize this dynamics? And what are those 
which can be implemented experimentally? Note that these questions need first to 
be answered for the classical atoms. We have to search for a configuation leading 
to complex classical dynamics. And experimentally, the dynamics in the classical 
limit must be characterized before considering the quantum system. In this paper, 
we address these questions, limiting our analysis to the classical limit. The quantum 
counterpart will be discussed in a future work. In section 2, we give some facts about 
cold atoms and optical lattices for those who are not familiar with this domain, and 
we discuss of possible lattices leading to complex dynanamics. Section 3 is devoted 
to the dynamics of atoms inside the wells, while section 4 deals with the dynamics 
of atoms traveling between several wells. Finally, we discuss in section 5 of the 
possible implementation of experimental measurements. 

2. Context: cold atoms and optical lattices. 

Cold atoms refer here to atoms cooled through a magneto-optical trap (MOT). 
The cooling is mainly obtained through an exchange of the momentum between an 
atom and a counter-propagating optical beam: while the absorption of a photon 
by the atom leads to a deceleration of the atom in the direction of the beam, 
the spontaneous re-emission of the photon arises in a random direction, and so 
does not change, in average, the atom velocity. To slow down atoms in 3D, three 
pairs of counterpropagating laser beams are necessary. Obviously, a moving atom 
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is decelerated by the photons traveling in a direction opposite to its own, but is 
accelerated by the photons traveling in the same direction as its own. But the 
frequency of these trap beams is detuned to the red of the atomic transition, so 
that, because of Doppler effect, the front photons are closer to resonance, and thus 
the deceleration process is more efficient than the acceleration one. This Doppler 
cooling process is coupled to an inhomogeneous magnetic field, which enhances the 
cooling process through the Zeeman levels splitting, and adds a restoring force to 
increase the atomic density of the cloud of cold atoms. MOTs lead typically, for 
Cesium atoms, to a 2 mm diameter cloud of 10 8 atoms at 5/xK. Such a cloud of cold 
atoms can exhibit spatio-temporal instabilities and chaos (Wilkowski et al. 2000, 
di Stefano et al. 2003, Hennequin 2004, di Stefano et al. 2004), but an adequate 
choice of the experimental parameters leads to a stable cloud, with atoms whose 
residual motion is the thermal agitation. 

A classical atom follows the motion equations of any classical object, and in 
particular the Newton's second law F — mr, where F is the force, r the position 
and m the mass (in the following, we take m = 1). When such an atom is dropped in 
a stationary wave, it undergoes a force F, the potential U of which is proportional 
to the wave intensity /, and inversely proportional to the detuning A between the 
wave frequency and the atomic transition frequency: 



Thus, atoms accumulate in bright (resp. dark) sites for A < (resp. A > 0). When 
the atoms are cooled with the MOT, the atomic density in these optical lattices 
is small enough to neglect the collisions between atoms, and so the only source of 
dissipation is the spontaneous emission. As spontaneous emission is proportional to 
//A 2 , it is relatively easy to build conservative optical lattices. Moreover, the clas- 
sical or quantum nature of the atoms in the lattice can be adjusted continuously, 
as it depends on the ratio between the temperature (or energy) of the atoms, and 
the depth of the lattice wells. For wells deep enough as compared to the atom tem- 
perature, the quantum properties of the atom, and in particular tunnelling, vanish, 
and thus atoms can be considered as classical (Greiner 2001). In the following, we 
always consider classical atoms, as discussed earlier. 

The atom dynamics in the lattice depends on the dimensionality of the lattice. 
For example, in a ID lattice, atoms have only two dynamical degrees of freedom, 
and thus even if the potential is not harmonic, the dynamics cannot be complex. 
It is needed to add at least a periodic forcing in such a lattice to observe chaos. On 
the contrary, a 2D lattice can exhibit chaos, without external forcing. 

But the atom dynamics also depends on the lattice geometry, and numerous 
lattice geometries can be obtained, as e.g. a vertical stack of ring traps (Courtade 
et al. 2006), five- fold symmetric lattice (Guidoni & Verkerk 1999) or even quasiperi- 
odic lattices (Guidoni et al. 1997). In this paper, we will focus on the case of two 
orthogonal stationary plane waves with the same polarization. The configuration 
of the laser beams is shown on Fig. la. The total field is £ = cos kx + e 1 ^ sin ky, 
where x and y are the two space coordinates, qb a phase, k — 2n/X the wave vector 



F = — 
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Figure 1. a) Layout of the laser beams, b) Spatial distribution of the intensity in the (X, Y) 
space. Black corresponds to the minimum value (zero intensity), while white corresponds 
to the maximum. The dotted square delimits the elementary mesh of the lattice, and the 
white crosses are the saddle points. 

and A the wavelength of the laser beam. The intensity can be written as: 

I = cos 2 kx + cos 2 ky + 2a cos kx cos ky (2-3) 

where a = cos <f>. With the adequate normalization, the potential is 

U± = ±1 (2.4) 

where the explicit sign is that of A. When a = 0, the coupling between x and y 
disappears, and the problem becomes separable. In all the other cases, the coupling 
between x and y could induce complex dynamics. It is easy to see that in these 
cases the elementary mesh of the potential is turned of ir/4 as compared to the 
(x, y) axes, and thus it is natural to introduce the following new coordinates: 



X = kx + ky (2.5) 

Y = ky-kx (2.6) 

The intensity and the potential can now be written: 

I = U+ = -U- = 1 + a(cosA + cosy) + cosAcosF (2.7) 



Before studying the potential, let us concentrate on the intensity. As an example, 
Fig. lb shows the spatial distribution of the intensity for a — 0.5. The elementary 
mesh is indicated through the dotted line. Assuming a > 0, the intensity I has 
an absolute maximum 2 (1 + a) at coordinates (n2n, to27t), where m and n are 
integers. It has also a relative maximum 2 (1 — a) in (jr + n.27r, it + m27r). Once 
again, we see that a = is a special case because the absolute and relative maxima 
have the same height. Note that a = 1 is another special case, where the intensity 
at the relative maximum vanishes and, thus, is equal to the minimum value. In this 
special case, we have black lines along X = tt + u2tt and Y = ir + n2ir. We will 
not consider these cases in the following. On the other hand, the intensity goes to 
zero in (tt + n2ir, m2ir) and (n2ir, tt + m2ir). Two neighboring zeros are separated 
by a saddle point where the intensity has the value I = 1 — a 2 . It is important to 
note that these saddle points are on the bissectors, connecting on a straight line 



Article submitted to Royal Society 



How to characterize the dynamics of cold atoms in non dissipative optical lattices ? 5 



an absolute maximum to a relative one and again to the next absolute maximum. 
On the contrary, the saddle points do not stand on the straight line that connect 
two neighboring zeros. This will induce a huge difference in the dynamics of atoms 
in the lattice obtained for red detunings (A < 0), where the atoms are attracted 
in high intensity regions and the one for blue detunings (A > 0), where the atoms 
are repelled from these same regions. The bissectors are clearly escape lines for the 
atoms when A < 0, while it is not the case for A > 0. 

Optical lattices appear to be an exciting tool to study the dynamics of a conser- 
vative complex system, but how to characterize this dynamics in the experiments? 
What are the experimentally accessible quantities? The typical size of a lattice mesh 
is A/2, i.e. 426 nm for Cesium. As the diameter of a cold atom cloud is typically 
2 mm, the 10 8 atoms are dropped in 22 10 6 sites for a 2D lattice, which lead to 5 
atoms/site. At these scales, it is clear that there is no way to isolate an atom, and 
thus no way to follow its trajectory. Moreover, to see an atom, we need light, and 
thus the measure introduces a dissipation and destroy the atomic state. A typical 
measure consists in illuminating the atoms with a laser flash, and recording the flu- 
orescence of the atoms through a camera. This destructive measure gives snapshots 
of the atom distribution in the space. We examine in the following if it is possible 
to extract informations about the atom dynamics from this type of measurement. 



3. Dynamics of atoms inside the wells 

Before we search for signatures of the dynamics in the experimental measurements, 
let us investigate in more details what are the relevant parameters and character- 
istics of the atom dynamics in a lattice. To illustrate this approach, let us consider 
again the two lattices introduced in Section 2. Although these two lattices differ 
only by the sign of their potential, they are deeply different. [/_ has its wells where 
the light intensity is maximum, while U+ has its wells where the intensity vanishes. 
Let us denote Et the value of the potential energy at the saddle point of the inten- 
sity. Atoms, the energy E of which is smaller than the threshold Et, are trapped 
into one site because they cannot climb up to the saddle point. On the contrary, 
atoms with E > Et can travel between sites, if they move in the good direction. 

Inside a trap site, the energy of the atom plays the role of a stochastic parameter. 
Indeed, for low energies, the atoms remain located close to the bottom of the well, 
and their dynamics can be approximated by an harmonic motion. As the energy 
increases, the potential becomes more and more anharmonic, the nonlincarities 
increase, and the dynamics can become more and more complex. To be able to 
compare the behavior of atoms in different potentials, we take in the following the 
origin of the energy at the bottom of the wells, and normalize the energy so that 
Et = 1. The potential energy then takes a different form for red and blue detunings: 

u + = ih ^ 

U. = 2(1 + V (3.2) 

Let us now examine in details the dynamics of the atoms in our two potentials. 
The most relevant way is to look at the evolution of the Poincare sections as a func- 
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tion of the energy. Our phase space is 4-dimensional, with directions (x, Y, X, YJ , 

but because of the energy conservation, the accessible space reduces to a 3D sur- 
face. We choose to consider Poincare section at Y = with increasing values, and 

thus, Poincare sections are in the 3D space (x, Y, Xj , and they lie on a 2D surface 

Sp, which has the shape looking like a semi-ellipsoid. To represent the Poincare 

sections we can project them on the (X, Y) plane or on the more usual ( X, X^j 

plane. The latter shows the Poincare sections viewed from the vertex of the semi- 
ellipsoid. However here, because of the stiff sides of Sp, the projection in this plane 
leads to a confuse map, as many curves are projected at the same location, and 
thus are superimposed. On the contrary, the projection on the (X, Y) plane gives 
more details, and thus in the following, we often choose it. However, let's keep in 
mind that we look at a lateral projection of a "bell" , and thus that we superimpose 
its front and rear faces. 

As pointed out before, because of the normalization we choose for the energy, 
the form of the potential energy differs in the cases of blue or red detuned lasers. 
We investigate each cases separately. 

In the case of red detuned lasers, the potential energy takes the form : 



9, 9/, (1 — cos X) (1 — cosF) 

U- = (1 - cosX) + cos Y)- i i-i— 2 -(3.3) 

(1 + a) 

witha;^ = (l + a) -1 (3.4) 

This potential appears to be the sum of two simple pendula coupled through the 
third term. The frequency for oscillations with small amplitude is the same for the 
two directions. This degeneracy together with the coupling term leads to a strong 
synchronization of the motion in the two directions (Bennet et al. 2002). However, 
in contrast with the Huygens clocks, we do not have any dissipation here, so the 
frequency locking occurs in a more subtle way (Hennequin & Verkerk 2010). 

However, it is interesting to identify the resonances of the system. A very simple 
approach is to restrict the problem to the first anharmonic terms, similar to the 
undamped Duffing oscillator. We then look for a periodic harmonic solution in the 
form X = Xq cos (uit) and Y = Yq cos (uit + tp), with u> close to u>q. We drop terms 
at other frequencies (i.e. 3w) and we end up with six families of solutions. The 
first two are the trivial ones : motion along the X or the Y directions (Yq = or 
Xq = 0). The other four are obtained for X — Y and for ip = 0, w, ±tt/2. For a 
given energy E, the relations giving Xq and ui are not simple, and it is beyond the 
aim of this article to write them explicitly. For the large amplitudes considered in 
the following, the motion is not any more harmonic and we cannot keep only the 
lower order terms, but the main result remains: we have six periodic trajectories, 
leading to points in the Poincare section (except for the trajectory Y = 0, that we 
cannot catch in a Poincare section at Y = 0). In the 3D space, these points have 
the coordinates (0,-Y ,0), (±X o ,-X ,0) and (0, -X^,±uiX^). 

In Fig. 2, we show the dynamics in the U— potential for different normalized 
energies in the case of a = 0.5. These results have been obtained through numerical 
resolution of the equations of motion which are derived from the potential (|3.3j) . 
without the addition of any random quantity. All the described behaviors are thus 
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Figure 2. (X, Y) (left) and (^X, X^j (right) Poincare sections of the atomic dynamics in 
the U- potential, (a) and (b): E = 0.80; (c) and (d): E = 0.88; (e) and (f): E = 1.00. 

deterministic. For each value of the energy, we project the Poincare section on the 
(X, Y) plane (left figures) and on the (x, Xj plane (right figures) . For low enough 

energies (e.g. E = 0.8 in Fig. 2a and 2b), we see four distinct domains separated by 
an X-shaped separatrix. In each of these domains, the Poincare section is cycling 
around one of the non-trivial resonances found above. As the motion along X and 
Y is governed by the same frequency, and because of the coupling between these 
two pendula, a synchronization between the two directions occurs, through a phase 
locking between the two motions. The corresponding behavior can be described 
as mainly a U) periodic cycle perturbed by small sidebands (Hennequin & Verkerk 
2010). 

The dynamics in U— evolves only slightly when E is increased. The Poincare 
surfaces are always organized around the separatrix delimiting 4 areas. In each area, 
the nature of the motion is the same, namely phase locking between the motions 
in the X and Y directions. Chaos appears close to the separatrix for E ~ 0.88 
(Fig. 2c and 2d), but it remains marginal, even when E = 1 (Fig. 2e and 2f). This 
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Figure 3. (X,Y) Poincare sections of the atomic dynamics in the U+ potential, (a) 
E = 0.4, (b) E = 0.63, (c) E = 0.80, (d) E = 0.88, (e) E = 0.93 and (f) E = 1.00. 

very small extent is due to the original degeneracy of the frequencies of the coupled 
pendula and to the strong coupling between them (Hennequin & Verkerk 2010). 

For blue detunings (A > 0), the bottom of the well corresponds to / = 0, i.e. 
(X = 0, Y = 7r) sites. For the sake of simplicity, we shift the origin in Y by it, in 
order to have a trapped motion centered at the origin. Thus, we can write : 



9 / 9 , \ (1 — cos X) (1 — cosF) 

U+ = u> 2 0X 1-cosX +^or 1-cosy (3.5 

(1 — a z ) 

with ujq X — (1 + a) 1 (3-6) 



= (I-")" 1 (3-7) 



U 0Y 



Once again, this potential appears to be the sum of two coupled pendula. But now, 
the two frequencies for oscillations with small amplitudes are different: for the value 
a = 0.5 chosen here, the ratio y/3 of these two frequencies is irrational. 

For very small energies (Fig. 3a), the dynamics consists essentially in a regular 
motion around the bottom of the well, along a quasiperiodic trajectory with fre- 
quencies lox and uy close to coqx and ujqy- At the top of Fig. 3a, Poincare sections 
are those of atoms, the motion of which is essentially along the X axis. In Y = 0, 
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the trajectory is a periodic cycle along the X direction (edge of the semi-ellipsoid). 
At the opposite, the periodic cycle at the bottom of the figure corresponds to the 
situation where the atomic motion is exclusively along the Y axis (vertex of the 
semi-ellipsoid). Note that the nature of the motion along these quasiperiodic cycles 
is deeply different from those described with A < 0. Indeed, as cj x and uj y are 
very different, no locking occurs. In particular, in the spectrum of the motion, the 
two main frequencies are close to loqx and ujqy- 

As the energy of the atom is increased, the atom can climb more and more in 
the well, the frequencies ujx and toy change because of the anharmonicity of the 
potential, but the dynamics does not change fundamentally until E ~ 0.6. At that 
point a new feature appears: a stable periodic trajectory shows up as a cycle close to 
the bottom of Fig. 3b, obtained for E = 0.63. In fact, for amplitudes large enough, 
the frequencies u>x and ujy depart so much from their initial values loqx and loqy 
that a new resonance appears at ujy = 2ujx- 

For higher energies, the uy = 2ujx resonance grows and comes closer to the 
centre of the figure and influence a non-negligible fraction of the trajectories. In 
Fig. 3c, for E = 0.8, the resonance is clearly visible in Y ~ —0.74. In the (x, Y, X^j 

space, its Poincare section consists in 2 points (superimposed in the projection of 
fig. 3c), explored alternatively by the trajectory. Around this point, the Poincare 
sections are a double closed loop. The corresponding quasiperiodic motion consists 
in a perturbed uy = 2ujx phase locked periodic cycle, where the perturbation 
consists in small sidebands of u>x and uiy in the spectrum. Thus the separatrix 
appears here to be the limit between this phase locked and the unlocked behaviors. 
The central domain and the two linked lateral domains (bottom left and right) 
correspond to the phase locking. The difference between these two domains is the 
relative phase on the motion along X and along Y. In the two other domains (top 
and bottom), there is no locking between the ujx and coy frequencies. 

In E = 0.8 (Fig. 3c), all the trajectories are still periodic cycles or quasiperiodic 
tori. When the energy is increased further, chaos appears at E ~ 0.88, starting in 
the vicinity of the separatrix (Fig. 3d). Then, it expands with some quasiperiodic 
islands remaining (Fig. 3e), but finally, for E = 1 (Fig. 3f), the only significant 
quasiperiodic domains are those around the X and Y periodic cycles. Around the 
locked periodic cycles, a narrow area with tori remains, but chaos appears really to 
be dominant. 

We have shown in this section that it is relatively easy to find two slightly dif- 
ferent lattices with fundamentally different dynamics. These two configurations are 
easy to reach experimentally, as they differ only by the sign of the detuning. It 
would be interesting now to examine how to measure experimentally these differ- 
ences, and if these differences have an impact on the dynamics of atoms when they 
jump between sites of the lattice. The next section deals with the latter. 



4. Dynamics of atoms visiting several wells 

To travel from site to site, an atom needs to have an energy E ^ 1, but it is not a 
sufficient condition. Only atoms with an adequate trajectory will effectively escape 
from a well. This implies that for a given energy E > 1, at least two classes of 
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Figure 4. (X, Y) plot of the trajectories of 100 atoms in the (a) U+ (E = 2.66) and (b) 
U- (E — 1.07) lattices. Each atom starts in the central mesh, and move during the time 
t = 10 6 , which corresponds to more than 10 5 periods of oscillation at the bottom of a well. 



atoms can exist: trapped atoms remaining in a single well, and traveling atoms, 
which escape the wells. In fact, the situation is more complex, as we will see now. 

Let us first examine the dynamics of traveling atoms in the blue case (A > 0). 
We are here interested by atoms with an energy 1 < E < 4. Indeed, atoms with 
E > 4 have an energy larger than the potential maximum, and thus they "fly" above 
the potential, and their trajectory is purely ballistic. On the contrary, the dynamics 
of the atoms with an intermediate energy consists in complex trajectories visiting 
a large number of sites, as in a random walk. As our model is fully deterministic, 
it involves in fact chaotic trajectories. Fig. 4a illustrates such a chaotic diffusion: 
it reports the trajectories followed by 100 atoms. Such a trajectory is in fact an 
alternation of oscillations inside wells and of jumps between wells. Here, we know 
that chaos dominates inside the wells, and thus the chaotic nature of the diffusion 
is not surprising. However, as we will see below, the existence of chaos inside the 
wells is not a necessary condition to observe a chaotic diffusion. 

To think of an experimental characterization of this chaotic diffusion, a simple 
way would be to characterize the diffusion function, and to evaluate a diffusion 
coefficient. Fig. 5a reports the distance covered by atoms of high energy (E = 2.66) 
as a function of time. With such high energies, all atoms travel between wells. They 
all follow a similar behavior, characterized by a diffusion over a distance of the order 
of 10 3 for the time interval of the figure. Although there is a small dependence of 
these curves as a function of the energy of the atom, the orders of magnitude remain 
the same for all energies 1 < E < 4. The only difference is that for lower energies, 
some atoms remain trapped in their well, and so a second group of curves appears 
with atoms remaining within a short distance (smaller than the mesh, i.e. 2tt) of 
their initial location. 

In the red detuned situation (A < 0), the maximum of the potential is at 
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Figure 5. Distance covered by 100 atoms as a function of time. Parameters are the same 
as in Fig. 4: in (a), U+ and E = 2.66; in (b), U- and E = 1.07. 

E = 1.33. As in the blue case, atoms with an energy E > 1.33 have ballistic 
trajectories, and atoms with 1 < E < 1.33 exhibit a diffusive chaotic behavior (Fig. 
4b). The origin of chaos now is clearly in the jumps between wells, as the dynamics 
in the wells is regular. And in fact, there is a main difference as compared to the 
blue case: the diffusion scale is larger by one order of magnitude, on the whole 
interval 1 < E < 1.33. We did not check if the slower diffusion originates effectively 
in the chaotic trajectories followed by the atoms inside the wells, but it would 
be interesting to check in a future work how these chaotic behavior could slow 
down the atoms. However, the difference of one order of magnitude in the diffusion 
speed reveals that the macroscopic behavior of atoms could effectively be used to 
characterize the nature of the dynamics in optical lattices 

But there is another important difference between the two lattices: in the red 
case a third regime exists, neither trapping neither diffusing. It is illustrated on 
Fig. 4b, where trajectories appear following the two bissectors. These trajectories 
correspond to atoms traveling along the escape lines of the lattice, as they were 
described in section 2. These atoms follow in fact a ballistic trajectory where they 
travel very rapidly along the bissectors. For example, in Fig. 4b, the ballistic tra- 
jectories reach 10 6 in all directions, while the diffusive atoms reach only 2 10 4 of 
the same units in the same time. Note that the ballistic trajectories we discuss here 
occur as soon as the threshold E = 1 is reach, and only along the escape lines of 
the potential. 

Fig. 5b shows the distance covered by the atoms as a function of time. We have 
now clearly three groups of trajectories: trapped trajectories at bottom, diffusive 
trajectories for distances of about 10 4 , and ballistic trajectories at the top, for 
distances larger than 10 5 . The main difference as compared to the A > case is 
the cohabitation of ballistic and diffusive trajectories, even just above threshold. 
This put in evidence three specific time scales of the dynamics of atoms with a 
given energy, associated respectively with the trapped, the chaotic diffusive and 
the ballistic trajectories. 

In this section, we examined the dynamics of atoms, whose energy is large 
enough to escape the potential wells, but remains smaller than the potential max- 
ima. We focused on atoms traveling between wells, and found a different behavior 
for our two lattices. For the red lattice, atoms can be classified following two types 
of dynamics: the diffusive atoms exhibit a chaotic dynamics carrying them off their 
initial location; the ballistic atoms move away rapidly from their initial location. 
These behaviors are associated with two different time scale. But is it sufficient to 
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Figure 6. Number of atoms versus time in (a) U+ potential and (b) U- potential. 



identify these different regimes in a real experiment? We have also shown that the 
dynamics of atoms in the blue lattice is quite different, both for the diffusive regime 
and the ballistic one: the time scale of the former is one order of magnitude smaller, 
while the latter simply does not exist. Can we use these properties to characterize 
and distinguish experimentally the two lattices? These questions are discussed in 
the next section. 



5. Macroscopic signatures of chaos 

Our aim is to characterize the dynamics of the cold atoms in the optical lattice. 
As we are concerned by conservative lattices, we cannot hope to "film" in real 
time the atoms in the lattice, as it would introduce dissipation. Thus we have 
to find other techniques. As the specificity of each lattice concerns the traveling 
atoms, an experimental measurement aiming at characterizing these lattices should 
characterize these traveling atoms. 

Experimentally, the lattice is finite. So the traveling atoms will reach the edge of 
the lattice, and finally leave the lattice. Therefore a simple measure of the lifetime of 
the atoms in the lattice give informations about the trapped and traveling atoms. 
However, as there are several types of traveling atoms, the simple measure of a 
lifetime is not sufficient, and the lifetime curve itself, in particular its shape, must 
be analyzed. Thus we will plot now the number of atoms in the lattice as a function 
of the time. The shape of the curve and the lifetime itself should give informations 
about the traveling atoms, while the baseline gives the percentage of trapped atoms. 

In the experiment, all the atoms have not the same energy, but on the contrary, 
they exhibit a distribution of energies linked to their temperature. Thus the results 
shown below have been obtained by using a sample of atoms with an appropriate 
distribution of energy. 

Fig. 6a shows the number of atoms in the blue lattice versus time. To simulate 
the finite size of the lattice, atoms are removed as soon as they reach a distance 
Dl = 1000. The curve exhibits a plateau at short times, followed by a exponential- 
like decreasing to an asymptote. The plateau corresponds to the time needed by 
the first atoms to reach the edge of the lattice (in the simulations, all the atoms are 
supposed to be initially at the center of the lattice) . The decreasing corresponds to 
the diffusing atoms escaping the lattice, and the asymptote to the number of atoms 
trapped in wells. This behavior does not depend on the lattice size Dl, except that 
the life time of atoms increases. In fact, the distance Dl = 1000, i.e. about 150 
lattice meshes or 70 um for a Cs trap, is smaller by one order of magnitude than a 
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typical experimental realization. However, a value of Dl = 10 4 leads, for the data 
presented in Fig. 4a and 5a, to an almost flat curve, because the time series are not 
long enough. To reach such a distance, one should increase the evolution time by 
two orders of magnitude. 

Fig. 6b shows the number of atoms in the red lattice versus time, for Dl = 10 4 . 
The shape of the curve is qualitatively different from that obtained for the blue 
lattice. At short times, a fast decreasing appears, corresponding to the loss of the 
ballistic atoms. At large times, not visible on the figure, an asymptote is reached, 
corresponding to the trapped atoms. The intermediate decreasing correspond to 
the loss of the diffusing atoms. Note that the decreasing appears to be more or less 
linear. In fact, the shape of this part of the curves is the sum of the diffusing losses 
of different classes of atoms differing by their energy. As a function of Dl , this sum 
can exhibit very different shapes, from an exponential-like shape, as in Fig. 6a, to 
an almost linear shape, as in Fig. 6b. 

Fig. 6 shows that the measure of the lifetime of atoms in a conservative optical 
lattice provides qualitative and quantitative informations about the nature of the 
lattice and the nature of the dynamics of the atoms in the lattice, in particular about 
the chaotic diffusion. Therefore the measure of the atom lifetime, in particular the 
existence of several characteristic times in the decreasing of the atom number, 
appears to be a signature of the chaotic dynamics of atoms in the lattice. 

6. Conclusion 

We have shown in this paper that optical lattices are a good toy model to study 
experimentally the dynamics of conservative systems, provided that relevant exper- 
imental measures are found to characterize this dynamics. In particular, we show 
that changing a simple experimental parameter can lead to two deeply different 
lattices, where atoms exhibit very different dynamical behaviors. We have shown 
that these differences exist both in the local dynamics of atoms inside a well, and in 
the non local dynamics of atoms traveling between wells. We searched numerically 
for signatures of these different dynamics in the experimentally accessible quanti- 
ties, and found that the measure of atom lifetimes in the lattice gives numerous 
informations about the existence and the type of chaotic diffusion of the atoms. 

It would be interesting now to characterize more precisely the diffusion function, 
as a function of the experimental parameters, in particular the atom temperature 
and the lattice size, and obviously to test these results on a real experiment. 

We considered in this paper only the dynamics of classical atoms, and it appears 
that this dynamics is more complex and more subtle than usually considered. Simple 
statistical analyses are not enough to fully characterize this dynamics, and more 
suitable tools are necessary. It is important now to think about the consequences 
of these results in the quantum regime, and in particular, to what could be the 
equivalent measures in the quantum world. 
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